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ABSTRACT 

We present i?-matrix calculations of photoabsorption and photoionization 
cross sections across the K-edge of Mg, Si, S, Ar, and Ca ions with more than 10 
electrons. The calculations include the effects of radiative and Auger damping by 
means of an optical potential. The wave functions are constructed from single- 
electron orbital bases obtained using a Thomas-Fermi-Dirac statistical model 
potential. Configuration interaction is considered among all states up to n = 
3. The damping processes affect the resonances converging to the K-thresholds 
causing them to display symmetric profiles of constant width that smear the 
otherwise sharp edge at the photoionization threshold. These data are important 
for modeling of features found in photoionized plasmas. 

Subject headings: atomic data - atomic processes - line formation - X-rays: 
spectroscopy 
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Introduction 



The current fleet of X-ray telescopes ( Chandra, XMM-Newton, Suzaku) have revealed 
complex spectra from astrophysical sources such as active galactic nuclei (AGNs) and X-ray 
binaries. Besides the ubiquitous Fe absorption features, many other elements from C to Ni 



have been observed in multip 
spectrum of NGC 3783 from 



e ionization stages ; see, for example, the 900 ksec Chandra 



Kaspi et al. 



( 120021 ) . In order to accurately model systems 



producing such spectra, atomic data is needed for all of these species, particularly the 
resonance structure leading to the absorption features. Due to the wavelengths of these 
features, the relevant atomic data involves the K-shell of species from C to Ni where the 
currently available data has been of limited accuracy. Over the last several years, our 
group has worked on providing this data. We have calculated energy levels, wavelengths, 
Einstein ^-coefficients, Auger rates, and phot oionization cross sec t ions. Complete data sets 



have been computed for all ion stages of 



Mendoza et al 



NJ 



2004 



Garcia et al 



Bautista et al 



2004; 



e (IBautista et al. 



Kallman et al 



2003 



'aimer i et al. 



20041 ) . O (IGarcia et 



al. 



20091 ) . There has also been recent work on the C 



( 120101 ). This work is an extension of 



Palmeri et al. 



( 12008h and 



sequence by 



2003a.b 



2005) , and 



Hasoglu et al. 



Witthoeft et al. 



(120091) 



where we present data for the Li-like to Ne-like ions of Ne, Mg, Si, S, Ar, and Ca. Here 
we continue these calculations for the 23 non-neutral, Na-like to Cl-like ions of the same 
elements. 

K-shell features of the ion stages covered here have been detected in X-ray spectra 
of astronomical objects, most often due to fluorescent line emission. Notable examples 
i nclude the detection o f Ka emission from the ions Si II-V in the X-ray binary Vela X-l 



( jPaerels fe Kahn 



20031 ) . This reveals the existence of a wide range of ionization states 
coexisting in this system, and can best be explained if the emitting gas is highly clumped. 
Ka lines from S IV- VI, Ar VI- VIII, and Ca VI- VII have been detected in the X-ray 
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binary 4U1700-37 (IBoroson et al.ll2003l ). The flux in these lines is correlated with the flux 
in the X-ray continuum, thus reinforcing the fluorescent origin. The K-edge of neutral 



Ne has been detected 



( jJuett fc Chakrabarty 



rom interstellar gas along the line of sight to bright X-ray binaries 



20031 ) and provides evidence for local enhancements in the interstellar 
Ne/O elemental abundance. The data we present here is important both for its direct 
effect on observed spectra via absorption features, and also through its effect on ionization 
balance calculations, which in turn affect many of the observables from these elements. In 
particular, it is the detailed resonance and edge structure, missing in previous calculations, 
which is most important. 

There have been a large number of calculations for these systems going back to the 
1920's, but they are confined to the non-resonant background and few cross the K-edge 
region. For the ions considered in this work, there is good agreement in the background 



cross section both a 



jove and below the K-edge and are well characterized by the fits of 



Verner fc Yakovlevl (119951 ). Measurements have been mainly confined to the neutral systems 
and for high photo n energi e s. A c omprehensive assessment and fitting of measurements has 



been performed by 



Veigelel jl973 ). 



There has been rece nt progress measur ing the res onance structure ne ar the K-edge 



of Li-like sys tems: C 3+ (IMiiller et al. 



Be-like C 2+ f lScullv et al. 



20091 ) and B 2+ (IMiiller et al. 



20101 ). as well as 



20051 ). R- matrix calculations performed as part of those works 



are in good agreement with the measurements apart from small discrepancies in some 
resonance positions and heights. With the advances of free electron lasers to produce X-rays 



( iKanter et al. 



20061 ). we expect it will be possible over the next several years to perform 



similar experiments for the species under investigation here. 

Tables of the total and partial cross sections accompany this arti cle as on-line tables 



Full data sets are also available through the xstar atomic database (IBautista &: Kallman 
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20011) 



Numerical methods 



The .R-matrix method e mployed here is the sa me as our previous K-shell calculations 
and is described in detail in Isautista et al.l ([2OO3L The R - matrix method is based on 



the close- coupling approxim a tion of 



following 



Burke et al. 



toll); 



Burke fc Seaton 



Berrington et al. 



(11974 . 



1978 



(1971) which is solved numerically 



19871 ). Due to the complexity 



of these M-shell ions, the scattering calculations are carried out in LS'-coupling except for 
the Na-like ions where we use the Breit-Pauli i?-matrix framework (bprm). The main 
limitation of the LS calculations is that the results are not fine-structure resolved. As we 
will show in the next section, the total cross section is not sensitive to the configuration 
(or coupling scheme) of the initial state. However, for modeling purposes, we would like 
to have cross sections resolved to the fine-structure level, so we split the cross section 

nj — i 

according to iRaul (|1976l) . Our wave functions were obtained using autostructure 



(Badnell 



1986 



( iRalchenko et al. 



1997) where the term (or level) energies have been adjusted to match NIST 



20081 1 when possible. For energies which are not available in NIST, most 



importantly th ose for the K-hole st ates, we adjust to match the relativistic Hartree-Fock 



calculations of 



Palmeri et al. 



(J2008|). 



As in our previous calculations, the effect of rad iative and spectator 



Gorczyca fc Badnelll ( 11996 



Auger 



decay on 



20001 ) . 



resonances is included using the optical potential of 
Here, the resonance energy with respect to the threshold acquires an imaginary component 
which is the sum of the radiative and Auger widths of the core. The radiativ e widths are 



computed within the i?-matrix calculations following 



Robicheaux et al. 



(119951 ). The Auger 



1 http:/ /heasarc. nasa.gov/lheasoft/xstar/xstar.html 
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widths are provided from the relativistic Hartree-Fock calculations of iPalmeri et al.l (120081 ) . 
These effects are referred to as radiative and Auger damping since they cause resonances 
in the photoionization cross section to become smaller in area. The damping also makes it 
necessary to calculate both the photoabsorption and photoionization cross sections since 
not every absorbed photon will result in a photoelectron. 



Berrington et al. 



(119951 ) for the inner region 



We use the i?-matrix computer package of 
calculation and the asymptotic region code STGBfOdamp (Gorcyzca & Badnell 1996; 
Badnell, unpublished) is used to compute the photoabsorption and photoionization cross 
sections including the effects of both radiation and Auger dampings. The target expansion 
for each ion includes the Is 2 2/ 8 3/^, Is 2 2/ 7 3l N+1 , and Is 2/ 8 3l N+1 configurations where 
iV+ 10 is the number of electrons in the photoionized (target) ion. Only configurations with 
at most a single 3d electron are included. The number of target states therefore can be as 
few as 43 levels (Na-like) or as many as 132 terms (P-like). We include enough continuum 
basis orbitals to span from threshold to at least three times the energy of the K-edge. This 
choice gives reliable cross sections up to at least 1.5 times the energy of the K-edge. Partial 
and total cross sections are calculated from all terms in the ground configuration of the 
parent ion. 

At this point, we would like to make a note about the accuracy of the K-edge energy. 
For neutral systems where measurements exist, we can be fairly confident of the edge 
position. However, for ions, we must rely on calculated positions. Due to the magnitude 
of these energies, differences on the order of a percent yield uncertainties in the edge of 
a few Rydberg. In the present calculations, the energy of the edge is determined by two 
quantities: the ionization potential cal culated by -R-matrix and the energy of the first 



K-hole target state which we take from 



Palmeri et al 



(2008). Since the ionization potential 



is known, we can (and do) adjust our cross sections by the difference with observations 
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(~ 0.2 Ry) when transferring our data to modeling codes. It is more difficult to judge the 
accuracy of the K-vacancy target state ener gies except to say tha t we expect the accuracy 



to improve as the charge state increases; see 



Palmeri et al. 



(120081 ) for an assessment of the 



accuracy. We estimate that the resonance and edge p ositions of the s ingly- ionized systems 



2010). We consider this to 



could be shifted by up to 1 Ry from their true values ( IGorczycal 
be the limiting error for these systems. The accuracy of the edge and resonance positions 
should improve significantly for the higher charge states. 



3. Results 

The cross sections reported here contain many of the same features shown in 
calculations from other ions and elements. Inclusion of radiation and Auger damping 
broadens resonances near the K-edge and makes them more symmetric. The damping is 
also important in separating the absorption and ionization cross sections. In Fig. (TJ we 
show these cross sections for the 4 S ground term of P-like Ca. There is little difference 
between absorption and ionization for the first, K/3 resonance which indicates that damping 
plays little role here. For the near-edge resonances, however, the photoionized resonances 
are nearly damped away entirely while the strongly-broadened absorption resonances lead 
to the characterstic smearing of the K-edge. Small resonances can be observed above 
the edge in both absorption and ionization. These features are either typical resonances 
converging onto a higher-energy edge or resonances belonging to doubly-excited states 
which are accessible via mixing. They are observed in all of the calculations, but rarely 
contribute much to the total cross section. 

In Fig. [21 we show the photoabsorption cross section for each of the Na-like ions. 
The behavior of the cross section across the iso-electronic sequence is quite regular. The 
background cross section above and below the K-edge increases with nuclear charge and the 
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width of the resonance region also increases with Z. However the backgro und cross section 



is kn o wn to be insensitive to the number of electrons in the valence shell ( IMissavage et al. 



1977 



Nasreen et al. 



19891 ). If we were to plot the cross sections for all iso-electronic 



sequences on this same figure, it would look much the same. The exception is below the 
L-edge for sequences with more than 12 electrons where the inner-shell (3s) electron has 
the same principal quantum number as the valence (3p) electron. In these cases, the 



insensitivity is lost (IPradhan et al. 



20091 ). 



In Fig. [31 we show the background c ross section at a fixed energy above the K-edge 



where results from 



Witthoeft et al 



(120091 ) are used for ions with 10 electrons or fewer. We 
find the background cross section for each sequence to be constant within 10%. Note that 
the choice of energy is arbitrary and is different for each sequence so no conclusions should 
be drawn about the relative cross sections between sequences. Nor does the choice of energy 
affect the agreement, we find better than 10% agreement between the cross sections at all 
energies where we can make comparisons. Only the Mg sequence shows a steady downward 
trend of the cross section with an increasing number of electrons, but the total drop is only 
10% so it is difficult to make any conclusions. The remaining sequences show no trend with 
number of electrons. Note that we are starting with the Li-like ions; if we also showed the 
cross section for photoionization of the He-like ions, we would expect to see differences since 
there is no inner-shell. 

Previous experiments and calculations for K-shell photoionization were more concerned 
with high-energy behavior than the detailed structure of the cross section near the K-edge. 
T hese results are generally i n good agreement and have been characterized by the fits 



of I Verner fc Yakovlevl (jl995l ). We want to compare our results with experiment, but 
measurements have only been taken for neutral systems not calculated here. However, 
we can take advantage of the finding above and compare our singly-ionized results with 
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the neutral measurements as a test of the background. In Fig. HI we show our cross 
sec tion for Cl-like Ar compar ed to the photoabsorption cross secti ons of neutral Ar given 



by 



Verner fc Yakovlevl (119951 ) and the fit to experiment al data by 



Vei geld ( 



Verner fc Yakovlevl (119951 ) at all 



1973). There 



is very good agreement between the present results and 
energies except at the K-edge where resonances give an enhancement to the i?-matrix 
results. The experimental data point at the K-edge seems to confirm this enhancement, 
although for the other seq uences the Veigele data at the K-edge are in good agreement with 
Verner fc Yakovlevl (119951 ). Except for the area in the immediate vicinity of the K-edge, 



there is similarly good agreement between all three results for the other elements. 



Our results are to be used in the xstar program ( jBautista fc Kallmanl 120011 ) for 
modeling photoionized plasmas. To prepare for this application, we convolve our data with 
a Gaussian and extrapolate to high energies. As these calculations are focused only on 
the resonance featu res near the K-edge, the cross section data from the Opacity Project 



( ISeaton et al 



1994 ) is used for the L-edge region. The two data sets are joined together by 



hand at an energy between the L- and K-edges. 

The energy-dependent width used for the convolution is AE/E = 10~ 3 which is 
representative of current detectors. The accuracy of i?-matrix calculations starts to degrade 
at high energies depending on the number of continuum basis functions included. When 
this happens, we start to see oscillations in the cross section. Therefore, the calculations are 
carried out to about 1.5 times the K-edge and then extrapolated to 1100 Ry (or 15 keV). 
The extrapolation assumes the asymptotic form, u{E) = AE~ P where the parameters, A 
and p, are determined by a least squares fit to the last ~100 points of the i?-matrix data. 
This method works well except for some weak partial cross section data where we start to 
see os cillations. To prevent unphysical extrapolation slopes, we enforce a lower limit of 8/3 



on p (I Cowan 



19811 ) . While the extrapolation is not as accurate for these weak partials, it is 
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not likely to affect modeling. 

Finally, after convolution and extrapolation, we remove unneccessary data points from 
each data set. Points which lie on a straight line between the two adjacent points within 
1% are removed. This test is repeated over the entire cross section until no more points are 
removed. Using this process, we can decrease the number of points in the cross section data 
from several thousand to a couple hundred, yet still maintain good accuracy with linear 
interpolation. In Fig. |5l we show the raw and convolved total cross sections for Mg-like Ar 
near the K-edge. The raw data contains over 11 000 points while the convolved data has 
less than 350 points. The convolved cross sections are not included with this paper, but are 
available by a request to MCW or xstar. The raw cross sections are available as electronic 
tables attached to this work. 



4. Summary and conclusions 



Total photoabsorption and photoionization cross sections have been computed for the 
K-shell of all non-neutral ions with 11 to 17 electrons for Mg, Si, S, Ar, and Ca. Partial 
photoionization cross sections have also been calculated for the same ions. Radiative and 
spectator Auger dampings are accounted for in detail and the energy region around the 
K-threshold was accurately treated for each ion. 



We find good ag reement in the back g round total cross section with experiment (jVeigele 



19731 ) and the fits of IVerner fc Yakovlevl (119951 ). The data provides term- resolved partials 



which, to our knowledge, have not previously been available. For the Na-like ions, the 
present results are fully level-resolved. 

All data are provided as on-line tables acc ompanying; this artic l e and can also be 



obtained through the xstar atomic database (IBautista fc Kallman 



20011) 1 and the 



Universal Atomic Database 



a 



The data set s provided her e 



rates reported in 



Palmer i et al- 



toge ther with the energy levels and radiative and Auger 



( 120081 ) will help modelers to carry out detail studies of K 



spectra of astrophysically abundant elements. 
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2 http:/ /heasarc. nasa.gov/uadb 
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Fig. 1. — Photoabsorption (thick, red) and photoionization (thin, blue) cross sections of 
P-like Ca near the K-edge. Color available in the on-line version of the text. 
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Fig. 3. — Background cross section at a fixed energy above the K-edge across iso-nuclear 
sequences. The photon energies for each sequence are: 128 Ry (Mg), 178 Ry (Si), 233 Ry 
(S), 300 Ry (Ar), and 375 Ry (Ca). Color available in the on-line version of the text. 
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Fig. 5. — Raw (thick, red) and convolved (thin, blue) cross sections for Mg-like Ar near 
the K-edge. The top plot shows photoabsorption and the bottom is photoionization. The 
convolution uses a width of AE/E = 0.001. Color available in the on-line version of the text. 



